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Abstract 

The single scalar field inflationary models that lead to scalar and tensor per- 
turbation spectra with amplitudes varying in direct proportion to one another 
are reconstructed by solving the Stewart-Lyth inverse problem to next-to- 
leading order in the slow-roll approximation. The potentials asymptote at 
high energies to an exponential form, corresponding to power law inflation, 
but diverge from this model at low energies, indicating that power law infla- 
tion is a repellor in this case. This feature implies that a fine-tuning of initial 
conditions is required if such models are to reproduce the observations. The 
required initial conditions might be set through the eternal inflation mecha- 
nism. If this is the case, it will imply that the spectral indices must be nearly 
constant, making the underlying model observationally indistinguishable from 
power law inflation. 
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I. INTRODUCTION 



A cosmological model with a low density of cold dark matter (^cdm ~ 0.3) and a 
dominant cosmological constant (^a ~ 0.7) is currently favored by recent analyses of cos- 
mological observations, in particular those of the cosmic microwave background (CMB) 
power spectrum (see for instance, f^-^) and high red-shift surveys of type la supernovae 
These analyses maximize the likelihood between the observed power spectrum and 
that of a given theoretical model by fitting values for various model parameters. This set 
of parameters includes quantities that characterize the initial conditions for the evolution of 
the density (scalar) and gravitational wave (tensor) perturbations. If these initial conditions 
are described by adiabatic, gaussian fluctuations, the observations strongly constrain the 
present-day universe to be very nearly spatially fiat 

The simplest casual mechanism for producing a spatially fiat universe, where large-scale 
structure originates from a primordial spectrum of adiabatic density perturbations, is given 
by the single field inflationary model 0. The early universe undergoes a rapid, accelerated 
expansion that is driven by the self-interaction potential energy of a real scalar field - 
the inflaton. We consider inflation driven by a single infiaton field throughout this paper. 
The field undergoes quantum fluctuations as it slowly rolls down its potential and these 
fluctuations generate a perturbation in the spatial curvature A primordial spectrum 

of tensor perturbations is also generated from the scalar field fluctuations [^|lOi. 

Until recently, the majority of CMB data analyses have neglected the possible effects 
of the primordial gravitational wave spectrum (For a recent review, see, e.g., Ref. @]). 
In the last few years, however, there has been a growing recognition that the role of the 
tensor perturbations deserves more attention when determining the best-fit values of the 
cosmological parameters [p|,|3|,p!T| . It has been found that if the tensor modes are included in 
the analysis, an extra degeneracy arises when fitting cosmological models to the CMB data 

Further motivation for including the tensor modes in the analysis arises from the possibil- 
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ity that a measurement of the CMB polarization may indirectly determine the gravitational 
wave contribution to the power spectrum [^,|T^. Such a contribution can be parametrized 
in terms of the quantity 



'^^a^' (1) 



representing the relative amplitudes of the tensor {At) and scalar {As) perturbations, where 
the constant, a, depends on the particular normalization of the spectral amplitudes that is 
chosen. 

Given the current accuracy of the observations, it is only possible at present to place 
a constant upper bound on the allowed range of r Looking to the future, the most 

optimistic expectation we have for constraining r is that the observations would favor a 



constant central value [TT|JT^. It is unlikely that a direct dependence of r on scale would 
ever be measurable. 

Motivated by these future expectations, therefore, it is important to determine the main 
features of the single field inflationary models that yield perturbation spectra resulting in a 
tensor-to-scalar ratio, r, that is precisely constant, at least over the range of scales accessible 
to observations. This is the purpose of the present paper. Of particular interest is the 
functional form of the corresponding infiaton potentials. Such a study is also relevant to the 
question of whether the infiaton potential can be reconstructed directly from observations 
and, consequently, provides important information about physics at very high energy scales 

HI- 

It follows immediately from the definition (|l]) that imposing the condition r = constant 
is equivalent to searching for models where the scale dependences of the two spectra are 
identical. In general, these scale dependences are parametrized in terms of the scalar {ns) 
and tensor {ut) spectral indices, respectively. These are defined in terms of the logarithmic 
derivatives of the corresponding normalized amplitudes: 

2 - dink ' 
ut _ dlnAr 

^ = T = ^hrr' 
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where k = aH is the comoving wavenumber when the mode first crosses the Hubble radius, 
dn = H~^, during inflation. 

We are therefore interested in the models that produce 

A(A;) = 5{k). (4) 

Since exact expressions relating the perturbation spectra to the inflaton potential are 
presently unknown, the standard approach is to expand the power spectra in terms of a 
set of parameters that describe the inflationary dynamics through the slow evolution of the 
Hubble radius or, in the case of single scalar fleld models, in terms of the slow rolling of the 
inflaton fleld along a nearly flat potential (see Refs. |l7|, p!^J18| and references therein). To 



leading order in these expansions, the inflationary models that have precisely constant ratio, 
r, are precisely the power law inflationary models driven by an exponential potential [[13 



For power law inflation the spectral indices are constant and equal to each other. However, 
as we discuss in the next section, to next-to-leading order, the indices can vary with scale 
whilst being equal to each other at each value of k. Speciflcally, we solve the next-to-leading 
order Stewart-Lyth inverse problem (SLIP) ||TB| under the condition that the spectral in- 
dices are equal at any scale. We flnd that for this case, power law behavior is a repellor 
rather than an attractor of the inflationary dynamics, in contrast to the lowest-order anal- 
ysis of Hoffman and Turner |]l6l. This difference arises because different approximations 



in the slow-roll analysis are considered. Nevertheless, if the underlying inflaton potentials 
derived at next-to-leading order are to produce successful inflation under the condition 
that r = constant (A = 5), a strong flne-tuning of the initial value for the inflaton fleld 
is required. The eternal inflation mechanism provides a natural way of obtaining the 



required initial values, although in this case, the scale dependence of the spectral indices is 
strongly suppressed and, consequently, the difference between these models and the power 
law model is effectively erased. 

The paper is organized as follows. We present the equations that determine the perturba- 
tion spectra in Section II. In Section HI, we proceed to solve these equations and reconstruct 



the corresponding inflationary potentials in parametric form. We conclude with a discussion 
in Section IV. 



II. PERTURBATION SPECTRA 

When determining the cosmological parameters the primordial perturbation spec- 

tra are often expanded as 

\nAl{k) = InAliK) + A{K) In A + In^ ^ + ■ ■ ■ , (5) 

A;^, 2 dink 

InAUk) = InAlik.) + In^ + lf^\k=k. In' ^ + ■ ■ ■ , (6) 

where k^: is a pivotal scale. The order where these expansions are truncated is determined 
by the precision of the CMB observations. Only the flrst two terms in each expansion are 
usually taken into account and the 'running' of the indices, dni/dlnk, is neglected. Even 
so, there are still four inflationary parameters that need to be fltted to the data: As{k^,), 
Axik^:), A{k^:) and 6{k^:). On the other hand, the total number of cosmological parameters 
in the analysis can be very large [^], and in view of the level of accuracy expected from 
forthcoming observations in the near future, it is possible that higher-order terms in the 
expansions and may have to be included as well. The large number of parameters 
that need to be considered allows degeneracies to arise when fltting a given theoretical model 
to the observational data. 

One way of reducing the number of degrees of freedom is to flt the ratio, r, instead of 
the tensor quantities directly. The relevant tensor modes can then be deduced from the 
definition of r, Eq. (|l]), and the next-to- leading order consistency relation for the class of 
single field models [0]: 



riT = -2 



--f-)'-(ns-l)(- 

a \aj \a 



(7) 



The bounds |^,|TT],|T2[ in the observed precision of the tensor contribution to the CMB power 
spectrum are such that even after a best-fit value for r{k^) has been deduced, the simplest 



approach to adopt when constraining the parameter space would be to assume that the 
value r{k^:) holds for all scales, i.e, r{k) = r{k^) = constant. As we have already seen, this 
corresponds formally to a model with A{k) = 6{k). Consequently, this procedure provides 
further motivation for analyzing the case where r is precisely constant, because this allows 
us in principle to gauge to what extent the accuracy of expansions and (|^) is sensitive 
to the value of fc*. 



In single-field inflation, the indices, A(fc) and S{k), satisfy the SLIP equations [|T5 



2Ceie\ - (2C + 3)eie\ - ei + + ei + A = , (8) 
2(C + l)eiei-e?-ei-<5 = 0, (9) 

to 'next-to-leading order' in the slow-roll approximation, where C = —0.7296 and a cir- 
cumflex accent denotes differentiation with respect to the variable r, defined such that 
dr = dlnH^. The first 'horizon flow function' is defined as IIT8|,|T5 



dlnc/n 3T 

where T = </)^/2 represents the kinetic energy of the inflaton field, and = ln(a/ai) is the 
number of ^foldings of inflationary expansion since some initial time, ti. (Note that the 
number of e-foldings is usually counted backwards in time. Here, we count it forward, i.e., 
N{ti) = 0). In general, Eq. (|lOD measures the logarithmic change of the Hubble distance per 
e-folding, or equivalently, the contribution of the inflaton field's kinetic energy relative to 
its total energy density. Inflation proceeds for ei < 1 (a > 0) and the weak energy condition 
for a spatially flat universe is satisfied for ei > 0. By combining Eqs. (H) and @ we deduce 
that 

6-A = 2Ceie\ - (ei + l)ei . (11) 
In the next Section, we proceed to find solutions to Eq. (0). 
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III. THE MODEL 



When determining the single field models that lead to r = constant we are interested in 
solving Eq. (|ll|) when Eq. (^) is imposed. Firstly, we remark that one solution to Eq. (p^D 
with this condition is that ei, A and 6 are all constants, corresponding to the power law 
infiationary model |]TB| . Moreover, power law infiation is the unique solution if second-order 
terms in Eq. (|TT]) are neglected. More generally, a first integration of Eq. (^) yields 



(12) 



when Eq. (§) is satisfied, where i? > is an integration constant. Substituting this result 
into Eq. (P) we find that 



Eq. ([13D is plotted in Fig. |l| for different values of B. 



(13) 




FIG. 1. The variation of the scalar and tensor spectral indices, 5 = A, as a function of the first 
horizon flow function ei for different values of the constant B. 



Now, it has been shown in Ref. |T^ that to next-to-leading order, the constraint 5(ei) < 
must be satisfied for any value of ei during infiation. However, as indicated in Fig. Eq. 
(p!3|) is a generalized inverted parabola and always has a positive maximum. This implies 
that 5(ei) > for some range of ei. Indeed, the positive maximum is located in the interval 



ei G (0,1/5). This would seem to indicate that the condition for the spectral indices 
to be equal (but not constant) might not be self-consistent at this level of the slow-roll 
approximation. However, the point is that 5 < is necessary only during inflation and, in 
effect, this condition further constrains the allowed range of values that ei may take for the 
analysis to be self-consistent. Indeed, we may specify the parameter, B, in such a way that 
the values S^e^rnax)) and e^max) are negligible, where 



^l(max) „ 



2exp(-^) 



(14) 



B{C + l) 

is the value of ei at the maximum value of 5(ei) and Cw is the Lambert W function [Q. As 
we shall see later, 1/ B determines the range of power law behavior of A. Thus, taking into 
account the range of ns consistent with analyses of CMB observations we can choose 

B = 100. In this way, ei(max) ~ 0.00025 can be identified as the lower critical value of ei 
where the precision in our calculations is consistent to next-to- leading order in the slow-roll 
approximation, with S{ei(max)) ~ 0.00009 then being consistently negligible. 
Before proceeding, let us first analyze the quantity 

since this term plays an important role in the forthcoming calculations. Substituting 6, as 
given by Eq. (pTS]), into Eq. (plSf ) implies that 

In order to obtain analytical results, it is necessary to approximate the denominator of Eq. 
(p!6|) in the different limits where ei ^ |ln(i?ei)| and ei ^ |ln(i?ei)|, respectively. In the 



former case, Eq. (0) can be be expanded as the series 

"ln(5ei) 



C ^ 1 ln(5ei 



ei 



(17) 



if ei is sufficiently different from zero. However, since ei is typically small during infla- 
tion, B must be appropriately tuned if this expansion is to be valid. Consequently, the 

8 



approximation is only appropriate over a very narrow range of e\. Indeed, a numerical in- 
vestigation confirms that the relevant range of ei is negligible for realistic values of B and 
so we regard it as unphysical. For example, \i B = 100, this approximation is valid only for 
ei G (0.009999,0.010001). 

The more interesting case is the limit, ei -C |ln (i?ei)|, where the following approximation 
for expression ( |T6|) applies: 

C \ 1 K ei 



Q 



1 



ei 



(18) 



+ ly eiln(5ei) 1 ln(5ei) 

Firstly, if ei < 1/-B, this implies that ei is very small and the absolute value of \\y{Be\) will 
be very large|^. Secondly, if ei > 1/-B, then ei G (l/S, 1] implies that ln(_Bei) > ei if S is 
sufficiently large. Note that the value of B specifies the central value of e\ during infiation 
and this is expected to be less than 0.1. With this in mind, a numerical test implies that 
the expansion is valid if ei > 0.112 for i? = 10 and e\ > 0.0101 if i? = 100. Furthermore, 
one may verify that this approximation is suitable even for ei close to 1/B. For example, 
when B = 100, the reliable intervals of ei are (0,0.0095] and [0.0101, 1], respectively. 

We now employ this approximation to find semi-analytical expressions for the spectral 
indices. The wavenumber, k = aH, at horizon crossing is evaluated as a function of ei by 
solving the first-order differential equation 

(C + l)(ei-l)ei-e2_ei-5 = 0, (19) 

where ei = dei/d\nk . Using approximation (|18D and integrating implies that 



ko 



+ 



_21n'(5ei) \n{Be 



B 



+ -^Si{l, t2 ln(5ei)) - ^S,{1, ^3 ln(5ei)) - In |ln(Sei) | 



(20) 



where 8i{n,x) represents the exponential integral [^], a =F corresponds to ei < 1/B and 
ei > 1/B, respectively, and a sub-dominant term has been consistently neglected. The 



^We are assuming implicitly that i? ^ 1 in this discussion. 
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integration constant is denoted hyf\ k^. 
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FIG. 2. The variation of the scalar and tensor spectral indices, A = (5, as a function of the 
Hubble radius crossing wavenumber, k. The upper branch corresponds to ei < 1/B and the lower 
branch to ei > 1/B. The power law solution, corresponding to ei = 1/B, is represented by the 
dashed line. 

The dependence of the spectral indices, 6{k) = A{k), on comoving wavenumber is pre- 
sentedQ in Fig. ^ It is observed that when the inflationary dynamics results in scalar and 
tensor perturbation spectra that satisfy r = constant, the power law inflationary model, 
as represented by the dashed line, is a repellor rather than an attractor to next-to-leading 
order. 

It is necessary to analyze the dynamics of the parameter, ei, in order to understand this 
behavior more fully. The solution to Eq. is given by 



r = exp <^ 2C 




B 



1 




+ 4^,(1, T21n(5ei)) 



253 
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(21) 



when the approximation ([T8|) is valid. Eq. ( pTD is plotted in Fig. 



^Integration constants are labeled by a subscript in the corresponding variable in what follows. 
^Hereafter the plots in this paper are consistently drawn using eK^^ax) &s the lower value for ei. 
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FIG. 3. The variation of ei as a function of r — tq = In H^. The flow of time along each branch 
is represented by the small arrows. 

It follows from this figure that the two branches in Fig. |^ correspond to initial values of ei 
that are greater than or less than 1/B, respectively. Since dr/dt < 0, the solutions diverge 
from the solution given by ei = 1/5 as cosmic time, t, increases and this is illustrated in 
Fig. ^by the direction of the small arrows near to each of the branches. In order to emphasize 
the (inverse) asymptotic behavior in the neighborhood of the unstable fixed point ei = 1/B, 
the initial values for ei for this figure were chosen to be extremely close to 1/B. However, it 
can be seen that any small deviation of the initial value from 1/5 is exponentially amplified 
and consequently, the time interval that the solution spends inside any given neighborhood 
of ei = 1/-B is exponentially suppressed. This implies that, unless the initial value of ei 
is extremely close to 1/B, this parameter will immediately move out of the [0,1) interval 
and consequently, this leads to a very rapid and undesirable end to the inflationary period. 
Observe that, in these cases, Tq determines the energy scale, Hq, at the end of the inflationary 
era. Returning to Fig. |^, such a fine tuning of the initial conditions for ei is translated into 
the requirement that a sufficiently large value of the integration constant, /cq, must be chosen 
in order to have near power law behavior in the appropriate range of scales. 

Further support for this conclusion regarding the necessary fine-tunning of the initial 
conditions is given by considering the variation of the e-foldings, N, with respect to ei. 
The lower limit, > 60, must be satisfied if inflation is to resolve the horizon and flatness 
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problems [0 and the dependence of the number of e-foldings on ei is determined from the 



differential equation, = — l/2ei ||2^ or, equivalently, by evaluating the integral 

dei 



N = -{C + l) f ,f\, +^o 
J ei + ei + 6 



(22) 



In our case, the integration of Eq. ( ]2^ ) yields 



N 



C 



1 



21n^(5ei) ln(5ei)J ' ln(5ei 



e? + — ^ + ln|ln(5ei) 



+ T In(Sei)) - t2 ln(5ei)) \ + , 



(23) 

and the graph corresponding to Eq. (|23D is shown in Fig. ^. As expected, the relative 
expansion rate of the universe is larger as the initial value of ei becomes closer to 1/B. 
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FIG. 4. Illustrating the variation of the number of e-foldings of inflationary expansion, A^, as 
a function of ei. Time flow along each branch is represented by the small arrows. 



The scalar field and its self-interaction potential are given in terms of ei by [|19 

1 



(3 - ei)exp [r(ei' 



2(C+ 1) 



eidei 



+ 



(24) 
(25) 



/2/t J el + ei + 6 

and the potential is shown in Fig. ^ where Vq = K~^exp(ro), k = Sir/mpi is the Einstein 
constant and mpi is the Planck mass. The right-hand branch of the potential has a minimum 



at e\"''''"' ' ~ 0.79, not shown in the figure in order to allow for the observation of details in 



(upper) 
-1 

a more realistic range of values for ei. 
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Unfortunately, even under the approximation ([T8|), it is not possible to analytically inte- 
grate Eq. (|25| ) for the infiaton field. Nevertheless, this integral can be performed numerically 
and the result is presented in Fig. 0. 




FIG. 5. The infiaton potential as function of ei. 







0.92 


0.94 


0.96 


0.98 


0,1 




-0.2- 


^"^^ 












-0.4- 
















t 


A 










-0.6 














t 


ir 












t 


T 










-0.8 
















t 


t 











FIG. 6. The infiaton field as function of ei. 

In Eq. (^), a sign for the square root of ei must be chosen and this specifies the signs of 
<p and dV/ d(j). We have assumed that > and dV/ d(f) < and so for consistency, one can 
only consider those intervals of ei where this is valid. With the above results it is possible 



to check that the necessary criteria [|l^] 



13 



are fulfilled only for Ci G (0, l/B) and ei G e\ 



{upper)-! 



(26) 



To conclude this section, therefore, 



the potential as function of the infiaton field is shown in Fig. |^. As can be observed in this 
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FIG. 7. The potential as function of the inflaton field. Vi denotes the branch where 
ei G (0,1/5), while Vn stands for ei G (1/5, e^"^^'^''^]. Larger values of the potential are as- 
sociated with ei being progressively nearer to the critical value l/B. In this limit, the model 
approaches the power law inflation behavior. 



figure, the required fine-tuning of ei is equivalent to a fine-tuning on the initial value of the 
infiaton field, i.e., a sufficiently high value for |0 — 0o| is required. In this plot, the branch 
where ei G (0, l/B) is denoted by V/, and Vu denotes the branch where ei G (1/5,6^"^^'^'^^]. 



IV. DISCUSSION 

Motivated by future prospects for measuring any possible contribution of the primordial 
tensor (gravitational wave) perturbations to the CMB power spectrum, we have investigated 
the class of infiationary models that result in scalar and tensor fiuctuation spectra with a 
constant ratio of amplitudes on all scales, r = constant. This condition is satisfied if the 
spectral indices are equal. To lowest-order in the slow-roll approximation, such a condition 
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of equality implies that the indices must also be independent of scale and this corresponds 
to the power law inflationary model. However, to next-to-leading order in the slow-roll 
approximation, the indices can be equal but may also exhibit a non-trivial dependence on 
comoving wavenumber. Under a self-consistent approximation, we have determined the 
functional form of this dependence and reconstructed parametric solutions for the infiaton 
potentials that produce such spectra. 

We find that there are two different possible potentials. Surprisingly, for the specific 
ansatz we consider, the power law inflationary model is a past, rather than a future, attractor, 
in the sense that at high energies (early times) both potentials converge to the exponential 
model, but move away from this special case at low energies. This provides a counter example 
to the generic, lowest-order analysis of Hoffman and Turner [|16| who find that power law 
inflation is a future attractor for the inflationary kinematics. This difference arises because 
the approach of Ref. [|1^] employs a strong version of the slow-roll approximation to analyze 



the constraints on the inflationary evolution. In particular, the effects of the inflaton's 
acceleration in its equation of motion and its kinetic term in the Friedmann equation were 
neglected. Moreover, the expressions relating the spectral indices and the ratio of the scalar 
and tensor amphtudes to the potential were truncated to lowest-order. In the present 
work, we have relaxed these restrictions and employed all of the available information on 
the infiaton dynamics and truncated the expressions for the spectral indices to next-to- 



leading order. Though the result of Ref. should contain some of the essential features 



of the dynamics, the highly non-linear nature of the next-to-leading order expressions leads 
to important deviations from these lowest-order results. Indeed, if the strong slow-roll 



approximation of Ref. [16 had been invoked in the present analysis, power law inflation 



would have emerged as the unique solution. However, the next-to-leading order dynamical 



analysis of Ref. indicates that in the reduced phase spaces for the evolution of ei, 
there exists a saddle point in the region where ei has interesting values and A < 0. This 
implies that with respect to cosmic time (recall that dr/dt < 0), attractor-like behavior will 
be characteristic only of those trajectories that are very close to the unstable separatrices. 
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Likewise, the saddle point acts as a repellor for those trajectories that are closer to the stable 
separatrices. Inflationary dynamics that does not correspond to a power law attractor was 
also found for the solution of the second case analyzed in Ref. . 

One consequence of the behavior described in the present work is that a strong fine- 
tuning of the initial value of the infiaton field is required if this model is to produce spectra 
consistent with observations. This is different to what typically arises in inflationary cosmol- 
ogy, where any sensitivity to initial conditions is washed out by the accelerated expansion. 
If we consider Figs. |l] and ^ once more, it can be seen that the upper branch of 6 grows and 
becomes positive, thus indicating a breakdown in the next-to-leading order analysis. For the 
lower branch, 6 begins to evolve extremely rapidly, probably indicating that the 'running' of 
the spectral index, dA/dlnk, becomes too large. Either way, observational constraints are 
difficult to satisfy. Thus, the potential in the region open to observation must be sufficiently 
close to the exponential (power law inflation) model and consequently the energy of the field 
initially stored in its potential must be sufficiently high. 

One way to satisfy this requirement is through the eternal inflation mechanism, where 
(large) quantum fluctuations in the infiaton can cause the field to diffuse up its potential 
P^ . In general, the condition for eternal inflation to arise is that [ P^ , P^ , 

^^>1, (27) 
vrm|,iei(i) 

for a given value of ei(tj). Recalling that dr = dhiH'^, we can rewrite this condition as 

r - To > In (7rei(i)mpi) (28) 

and consequently, t — tq must be sufficiently large for eternal inflation to proceed. On 
the other hand, it is clear from Fig. |^ that this is equivalent to requiring that the initial 
value of ei be close enough to 1/5. In this case, however, the model would effectively be 
indistinguishable from that of the power law model. Moreover, it could be argued that this 
restriction on ei{ti) itself represents a fine-tuning. Nevertheless, the attractive feature of 
eternal inflation is that only a small region of the universe need satisfy Eq. (pSj) for the 
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process of self-reproduction to start and continue indefinitely. If, as is commonly assumed, 
ei{ti) is randomly distributed, then at any given time in an inflationary universe described by 
this model, there is a finite probability for the existence of a region satisfying the necessary 
condition. 

In conclusion, therefore, our analysis indicates that it is difficult, from a theoretical point 
of view, to obtain spectra where r is truly constant if the spectral indices have a non-trivial 
scale-dependence. This particularly applies in the case where the slow-roll condition is not 
necessarily satisfied. Moreover, the potential must be close to the exponential form over the 
range of infiaton values accessible to observations. This implies that the observations would 
not be able to discriminate between the models discussed above and power law infiation. 
Consequently, if a running of the scalar spectral index is eventually favored by future ob- 
servations, it is likely that the underlying theory would be much more complicated and, in 
particular, would also result in a non-trivial scale-dependence for the ratio of amplitudes. 
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